Bunching and anti-bunching of localised particles in disordered 
media 



O 



C/3 



I 

o 



> 



O 



Frank Schlawin, Nicolas Cherroret and Andreas Buchleitner 

Physikalisches Institut, Albert-Ludwigs-Universitdt Freiburg, Hermann Herder Strafie 3, D-79104 Freiburg, Germany 



PACS 42 . 25 . Dd - Wave propagation in random media 
PACS 42 . 50 . -p - Quantum Optics 
PACS 72 . 15 . Rn - Localisation effects 

Abstract - We consider pairs of non-interacting quantum particles transmitted through a disor- 
dered medium, with emphasis on the role of their quantum statistics. It is shown that particle- 
number correlations measured in transmission are strikingly sensitive to the quantum nature of 
the particles when they undergo Anderson localisation, due to bosonic bunching and fermionic 
anti-bunching in the scattering channels of the medium. The case of distinguishable particles is 
also discussed. 



Wave transport in disordered media has been a sub- 
ject of intense research efforts for over 50 years. Interfer- 
ence phenomena such as Anderson localisation [1], coher- 
ent backscattering or universal conductance fluctu- 
ations [His] represent striking examples of the impact of 
disorder on wave propagation. While these phenomena 
can be explained on the level of wave mechanics, the im- 
pact of the quantised nature of the propagating wave or 
of the scattering medium has recently moved into focus of 
theoretical and experimental investigations. For instance, 
for classical light scattered resonantly off "quantum" dis- 
ordered media consisting of cold atoms, the atomic inter- 
nal structure has dramatic effects on the visibility of the 
coherent backscattering effect [TUH] . Conversely, the quan- 
tum nature of the propagating object can also come into 
play and lead to original phenomena, like in the context of 
cold atomic gases scattered off random optical potentials, 
where extensions of standard mean-field (wave-like) de- 
scriptions of matter waves allowed for instance to demon- 
strate Anderson localisation of the quantum excitations of 
Bose-Einstein condensates [B]. 

In optics, quantum statistics of multiparticlc interfer- 
ence effects move now into reach of multiple scattering 
experiments. Recent theoretical works notably include in- 
vestigations of photon noise and photon correlations in 
conservative [QfTO] , absorbing [TT]fT2| or amplifying [l2llT3] 
disordered media, or even of quantum entanglement in 
the presence of disorder [TIHTB] . Most of these previous 
works considered transport of non-classical light in diffu- 
sive media, i.e. where interferences are smeared out by a 



disorder average, but recently also the scattering of pairs 
of non-interacting quantum particles in one-dimensional 
disordered lattices was investigated under conditions of 
Anderson localisation [T7|. In particular, in this work the 
role of the quantum statistics of individual particles was 
pointed out, and bosonic bunching as well as fermionic 
anti-bunching effects were observed in particle-number 
correlations (known as Hanbury Brown-Twiss correlations 
in the absence of disorder [T5]). 

In the present Letter we present an original transport 
study of pairs of non-interacting quantum particles in 
waveguide-shaped media with three-dimensional disorder, 
where a cross-over from diffusion to Anderson localisation 
is present. Our purpose is here to clarify how the quan- 
tum statistics of the particles (bosons, fermions or distin- 
guishable particles) manifests itself in the particle-number 
correlation as the regime of transport is changed inside the 
medium. From both a numerical and an analytical study 
based on the random matrix theory of quantum transport, 
we find that, as long as propagation is diffusive in the dis- 
ordered medium, the correlation is hardly sensitive to the 
quantum statistics of the particles. In strong contrast, 
when Anderson localisation comes into play, the correla- 
tion exhibits completely different behaviour for bosons, 
fermions and distinguishable particles, due to bunching 
and anti-bunching effects in the scattering channels sup- 
ported by the medium. 

Model. — We consider a source producing pairs of 
monochromatic particles with identical wave number kQ. 
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These pairs propagate through a medium containing het- 
erogeneities randomly arranged in space and characterised 
by a scattering mean free path £, see fig. [TJ The medium 
has length L and cross-section A, and is assumed to be 
connected to perfect leads (not shown in fig. [T]). In such 
a geometry the incoming and outgoing wave fronts can 
be decomposed on a discrete set of ~ kg A transverse 
modes, each of them corresponding to a value of the trans- 
verse component of the wave vector compatible with the 
boundary conditions. The statistical properties of trans- 
port observables are independent of the transverse shape 
of the medium in the limit A'^ ^ 1 that we assume from 
here on. Particles' coincidences between two modes k ^ k' 
are detected in transmission. Such a single-mode detec- 
tion setup requires a detector area smaller than the typical 
size /N of a speckle spot, with R the distance between 
sample and detector. To quantify the coincident count 
statistics, we quantise the field before and after the disor- 
dered region, using the basis of the transverse modes [12] , 
and define the normalised correlation function 



C = 



■.UkTlk' 



{Uk) X iny 



(1) 



where fik denotes the particle-number operator in the out- 
going mode fc, (. . .) the quantum expectation value, : : the 
normal ordering of operators, and the overbar the average 
over a statistical ensemble of disordered samples. Since 
we are here primarily interested in the quantum statistics 
of the particles (fermions, bosons or distinguishable parti- 
cles) , we remain general in the following and do not specify 
how the two-particle state is initially distributed over the 
incoming modes. In the spirit of [M], we thus write the 
initial density matrix in the general form 



ai>aji 



(2) 



where the operator aj creates a particle in the transverse 
mode labeled by index z, and all sums run from 1 to N. 
Creation and annihilation operators fulfill the commuta- 



tion relations d^a 
Sij ± OjCii, with 



±a]d|, djdj 



±djdi and d]dj 



for bosons and — for fermions. Due 
to these relations, the tensor be chosen (anti-) 

symmetric under the exchange of the first and last two in- 
dices for (fermions) bosons. The prefactor 1/2 in cq. ([2]) 
guarantees the state normalization, tip = Xii j ''^ijij ~ 1- 
From the initial density matrix p, the mean particle num- 
ber in some outgoing mode k is defined as 



(fik) ^ tT[pcl{a')ck{a)], 



(3) 

where Cfc(a) = J^i^kiCii and cl{a^) = J^i^lt'^l ^^e the 
input-output relations connecting the annihilation (cre- 
ation) operator (c[.) of the outgoing mode k to the set 
of operators di (a|) of the incoming modes i through the 
(random) transmission matrix t of the disordered medium. 
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Fig. 1: A source produces pairs of particles characterised by a 
two-particle density matrix p. These pairs are scattered from 
a disordered medium with length L, mean free path £, and 
supporting A*' transverse modes. The coincidence rate between 
two outgoing modes k and k' is analysed in transmission. 



Substituting these relations into eq. ^ , we readily obtain 



E 



(4) 



Similarly, from the definition 

{■.hkfik'-) =^^[pc\{g))c\,{^)ck{a)ck'{a)] (5) 
of the particle-number correlation, we obtain 



f* 

k'j'^ki' 



tk'j'- 



(6) 



Note that at this stage eqs. Q and ([6]) hold for fermions 
and bosons, provided the tensor Wiji'j' is properly (anti-) 
symmetrised. 

The next step in the derivation of the correlation func- 
tion ([T|) consists in performing the ensemble averages. In 
eqs. (U) and (]§]), statistical properties of the disordered 
medium are encoded in the transmission matrix elements 
tki- In the limit N£/L = g ^ 1, which corresponds to 
a regime where particles are multiply scattered according 
to a diffusion process, these elements are normally dis- 
tributed over the unitary group jl9j . This property was 
notably used in |l3] to calculate the statistical distribution 
of {-.hkfik'') for photons. 

Here however, we aim at studying the behaviour of C 
for arbitrary values of the parameter g. For this purpose 
we make use of the unitarity and symmetry of the trans- 
mission matrix, which allows us to employ the polar de- 
composition [TOll^ 



tki 



Ukay^Va 



(7) 



where u and v are random unitary matrices and ti , . . . , r^r 
are the eigenvalues of t. Eq. (O is at the basis of the 
macroscopic description of disordered media fl9| . Physi- 
cally, the matrix v first distributes a field amplitude in- 
cident in mode i among the scattering channels of the 
medium, each having its transmission coefficient Tq (be- 
tween and 1). The matrix u finally recombines the trans- 
mitted field amplitudes in the outgoing mode k. 

For a medium with length much larger than its width, 
i.e in the limit L 3> \/A, ensemble averaging can be carried 
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out in two steps: first over the matrices u and w, which 
are uniformly distributed over the unitary group (isotropy 
assumption) , and second over the eigenvalues ri , . . . , t^v , 
whose joint distribution follows the so-called DMPK equa- 
tion [20]. Substituting the decomposition ([7]) into eqs. dH) 
and ([S]), and performing the average over unitary matrices, 
we obtain 



C 



1 Eg.fc Tan ± J2a 

2 iEa^f 



(8) 



where the positive sign refers to bosons and the negative 
sign to fermions. Eq. ([5]) is the first important result of 
this Letter. It shows that, when L ^ ^/A, the correla- 
tion function C is independent of how the incident state 
is distributed over the incident modes of the disordered 
medium, and only depends on the statistical correlations 
between scattering channels |14j . This property originates 
from the fact that when deriving eq. (|8]), we assumed each 
field amplitude incoming in a given mode to be equally 
distributed among the scattering channels of the medium 
(isotropy assumption). Another important comment con- 
cerns the prcfactor 1/2 in eq. ([5]), also found in [15[[TCj . 
For the particular case of a two-photon Fock state with 
Fano factor F = 0, this prefactor was shown to express 
itself as 1 + {F -l)/2 [TO]. Its value 1/2 < 1 therefore 
signals the sub-Poissonian statistics, here of an arbitrary 
incident two-particle state. 

We now comment on the structure of the mean corre- 
lation [numerator of eq. ([8])]. The latter consists of two 
terms. The first one, J2a b '''aTb, corresponds to the sum of 
all possible correlations between scattering channels. The 
second one, ± J2a ; corresponds to the sum of all possi- 
ble autocorrelations of scattering channels and is affected 
by the specific statistics of quantum particles: fermions 
cannot both propagate in the same channel, therefore this 
sum appears with a minus sign, and cancels the autocor- 
relation contributions in the first term (anti-bunching ef- 
fect). On the other hand, bosons tend to bunch in the scat- 
tering channels, as manifested by the positive sign in front 
of the second term. These bunching and anti-bunching ef- 
fects have important consequences on the magnitude of 
the correlation function C, as we now show. 

Numerics for bosons and fermions. Eq. is 

valid for any value of the parameter g = Ni/ L, provided 
L S> a/a. In order to clarify how bunching and anti- 
bunching effects would show up concretely in an experi- 
ment accessing the correlation function ([T]), we evaluated 
C numerically as a function of L/i, for a fixed number of 
modes, N = 10. For this purpose we computed each of 
the averages J^a "^a > J2a.b ^^^^ DMPK 
equation, by means of a Markov chain Monte Carlo ap- 
proach |22j . Such a numerical procedure was already used 
in the context of the conductance distribution of metal- 
lic conductors [53] and of intensity correlations of classical 
light in disordered media [24]. At fixed the study of 
C as a function oi Ljl allows us to probe all possible 
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Fig. 2: Correlation function ((8} for bosons (dots, blue online) 
and fermions (squares, red online), plotted as functions of Ljt, 
for A*' = 10. Diamonds (orange online) are the contribution 
Y.a'^IY.a^ to (7''°'°"'. Dashed horizontal lines are the an- 
alytical prediction ^ , and dashed curves are given by eqs. 
and (|12|l . Vertical dotted lines indicate the crossovers L ~ £ 
and L ~ N£ from ballistic to diffusion and from diffusion to 
localisation, respectively. 



regimes of transport in the disorder, from quasi-ballistic 
{L < £) to diffusive propagation {£ L <^ N£) and even- 
tually Anderson localisation {L ^ Ni). The results can be 
seen in fig. [5] for bosons (dots, blue online) and fermions 
(squares, red online). In the quasi-ballistic and diffusive 
regimes, C is very close to 1/2 for bosons and fermions, i.e. 
the correlation is hardly sensitive to the quantum statis- 
tics of the particles (the small difference observed in fig. 
[5] is discussed in the next section). This picture dramat- 
ically changes as L increases, at the onset of localisation 
[L ~ Ni): the bosonic correlation function starts grow- 
ing, while the fermionic counterpart decreases. Far in the 
localisation regime {L ^ Ni), C ^ 1 for bosons and 
C <C 1 for fermions. These radically different behaviours 
can be interpreted in the framework of the theory of ac- 
tive transmission channels |21[|25| by rewriting eq. ([5]) 

as ^ (1/2) X (Ea^,T^ + 2E.^)/(Ea^)' 

and C'"— = (1/2) X (Ea^b^) / (Ea^f ■ As the 
length L increases, more and more scattering channels be- 
come closed, i.e. their transmission coefficient becomes ex- 
ponentially small. Eventually, in the localisation regime, 
only one channel oq retains a fairly large transmission co- 
efficient. In this limit, we thus have (Ea^)^ — ''ao^: 

J2a '^a - 'Tao ^^'^ Ea^b"^^ - "W^j where Bq denotes a 
closed channel. Since the correlation between a closed and 
the open channel is typically very small compared to Tag , 
we have C'fCTmions ^ ^ other hand, we have also 

~ such that C^°^°^^ ~ (1/2) X IT/tZ:^. This 



ratio is very large because transmission fluctuations well 
exceed the mean transmission in the localisation regime 
[TO] . To get a better picture of the weight of each term ap- 
pearing in (^bosons (^fermions ^ ^^le different regimes 
of transport, we also show in fig. [5] the contribution of 
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channel autocorrelations to C^°^°^^, J2a'''a / J2a^^ ("^i^" 
monds, orange online). 

It is instructive to compare the results in fig. [2] for 
bosons (dots) with those obtained by Ott et ai, fig. 4 
of Ref. [TS], for the particular case of a two-photon Fock 
state. Our correlation function is insensitive to the form 
of the incident state, unlike in |15l , where the correlation 
function is constant for two photons incident in the same 
mode whereas it increases and saturates for two photons 
incident in two different modes. This difference relies on 
the different quantity appearing in the denominator of the 
correlation function of [TS] , (n^) (nj.) , instead of (uk) x (n'^,) 
in our case, eq. ([J). 

Analytical results. — In order to support our nu- 
merical observations, we now calculate explicitly C in the 
quasi-ballistic, diffusive and localisation regimes. First, 
when L < i, we have ~ 1 Va, such that 



I N^±N _ 1 ( 1 



(9) 



Eq. © is shown in fig. [2] for bosons (-I-) and fcrmions (— ) 
as dashed horizontal lines. The 1/A^ difference between 
the two types of particles is visible in the figure since we 
considered a finite value of N for the numerics. In practice 
however, ^ 1 and this difference is negligible. 

In the diffusive regime £ <^ L Ni, averages over 
transmission eigenvalues can be evaluated by making use 
of the method of moments introduced in [26] . The result is 
a perturbation expansion of C in powers of the parameter 
1/g = L/{Nl) « I: 



(10) 



In the regime of diffusive transport, bosonic bunching (-f-) 
and fermionic anti- bunching (— ) in the scattering channels 
are thus visible, but remain small as long as g 3> 1. 

Finally, far in the localisation regime L ^ N£, C can be 
calculated from the probability distribution of the trans- 
mission eigenvalues [12]. A straightforward calculation 
then yields: 



^bosons 




(11) 

(12) 



where we have explicitly introduced the localisation length 
^ = Ni. Eqs. ([TT]) and ((T2]) are consistent with the discus- 
sion above, namely C^"""'""" < 1 < C'''°"°"". Analytical 
predictions ([TT]) and (fT2]) are shown in fig. [2] as dashed 
curves, and agree well with our Monte Carlo simulations. 

Distinguishable particles. — To complete our dis- 
cussion, we finally consider the scattering problem of a pair 
of distinguishable particles. For this purpose, we equip 



them with an additional spin degree of freedom, indicated 
by the symbols f and |. With this strategy, we consider 
an incident state of the form 



P = 



E 



(13) 



In comparison with eq. note the missing prefactor 

1/2. This originates from the new commutation relation 
involving the two particles which have now different spins, 
^'^t^l^ ~ '^fe4,^'^t- Normalization of the state (jl3p again 
reads trp = J^ij = 1- 

The formalism developed above for indistinguishable 
particles remains essentially the same for distinguishable 
particles, but care has to be taken in defining the particle- 
number and correlation operators. Indeed, since one de- 
tects particles on output of the medium irrespective of 
their spin, these operators must now be symmetrised with 
respect to the two particles, i.e 



(14) 



for the correlation operator in the outgoing mode fc, and 
■.nkhk':= cl^cl,^Ck^Ck'i + cl^cl,^CkiCk'^. (15) 

for the coincidence rate operator between the outgoing 
modes k and k' . The input-output relations now read 



Cfct(i) = Y.^tkt'^^Ui) and c):.^(_^) = ^fc^'^'aJtCi) 
the spin labeling also in the transmission matrix elements). 
Quantum expectation values for fik and : fikfik' ■ are now 
given by 



and 



*t(i);.t 



[nk 



) = E W^3^'A^^' + E ^y ^' ^t/kf ' (16) 



iji' 



ijj' 



and 



■ fikfik'-) 



iji'i' 



[tttiU 



"^J* 3 \ ki ki' k' j k'j 



'■k'i''k'i'''kj''kj' 



(17) 

Note that in eq. (|16l) . the first (second) term on the right- 
hand side is nothing but the contribution of the particle 
with spin up (down) to the detected mean particle number 
in mode k. Similarly, the two terms on the right-hand side 
of eq. ()17p correspond to the two possible coincidence 
events, i.e. particle with spin up detected in mode k and 
particle with spin down detected in mode k' for the first 
term, and vice versa for the second one. Finally, we again 



decompose transmission coefficients according to t 



matrices u^'^ and 



''fca 



„ta) 



t(i) 

ki 



and carry out the averages over the 
This yields: 



^disting 



1 Ea.feTgTb 



(18) 



Comparing with eq. (jS]), we see that the bunching/ anti- 
bunching term has disappeared for distinguishable par- 
ticles, as expected. We computed C'^'*'''"^ numerically. 
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Fig. 3: Correlation function for bosons (dots, blue online) and 
distinguishable particles (triangles, green online), plotted as 
functions of L/£ for A'^ — 10. The upper and lower dashed 
curves are the analytical predictions (jlip and (|18p respectively, 
and the vertical dotted line indicates the crossover L ~ N£ 
from diffusion to localisation. 

using the same Monte Carlo approach as above to carry 
out the remaining averages over transmission eigenvalues. 
The results are shown in fig. [3] as a function oi L/£ (trian- 
gles, green online). For comparison we also show C't.osons 
(dots, blue online). As for bosons, the correlation func- 
tion of distinguishable particles increases with L/£. This 
is expected since at the onset of localisation, one scat- 
tering channel ao dominates transport and X^a b '''a.'Tb — 
T^^ ^ i^a^)"^ — Again, simple analytical results 

can be obtained in the quasi-ballistic regime of transport, 
where C''i"5''"g — 1/2, and in the diffusion regime, where 
C-disting ^ (^^^2) X [1 + 2/(15g2) + 0(1/^3)] . On the other 
hand, far in the localisation regime we find 



(19) 



This result is shown in fig. [3] (lower dashed curve), and 
agrees well with the numerical calculations. It also con- 
firms that the bunching effect of bosons is only fully visible 
in the localisation regime where 



C 



bosons 



(^disting 



(20) 



as also demonstrated by the numerical points in fig. |3l 

Conclusion. — We have shown that particle-number 
correlations measured on output of a disordered medium 
are extremely sensitive to the quantum statistics of lo- 
calised particles. This phenomenon reflects the bosonic 
bunching and fcrmionic anti-bunching in the scattering 
channels of the medium. An interesting extension of this 
work would be to analyse the scattering problem of many- 
particle quantum states. 
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